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Abstract 


Second-generation biofuels and bio-based products derived from lignocellulosic biomass are likely to replace 
current fuels derived from simple sugars and starch because of greater yield potential and less competition with 
food production. Besides the high aboveground biomass production, these bioenergy grasses also exhibit exten- 
sive root systems. The decomposition of root biomass greatly influences nutrient cycling and microbial activity 
and subsequent accumulation of carbon (C) in the soil. The objective of this research was thus to characterize 
root morphological and chemical differences in six perennial grass species in order to better understand root 
decomposition and belowground C cycling of these bioenergy cropping systems. Giant reed (Arundo donax), ele- 
phantgrass (Pennisetum purpureum), energycane (Saccharum spp.), sugarcane (Saccharum spp.), sweetcane (Saccha- 
rum arundinaceum), and giant miscanthus (Miscanthus x giganteus) were established in Fall 2008 in research plots 
near Gainesville, Florida. Root decomposition rates were measured in situ from root decomposition bags over 
12 months along with initial and final root tissue composition. Root potential decomposition rate constant (K) 
was higher in elephantgrass (3.64 g kg~' day~') and sweetcane (2.77 g kg~' day~') than in sugarcane (1.62 g 
kg" day~') and energycane (1.48 g kg! day '). Notably, K was positively related to initial root tissue total C 
(Total C), total fiber glucose (TFG), total fiber xylose (TFX), and total fiber carbohydrate (TFC) concentrations, 
but negatively related to total fiber arabinose (TFA) and lignin (TL) concentrations and specific root volume 
(SRV). Among the six species, elephantgrass exhibited root traits most favorable for fast decomposition: high 
TFG, high TFX, high TFC, high specific root length (SRL), and a low SRV, whereas giant reed, sugarcane, and 
energycane exhibited slow decomposition rates and the corresponding root traits. Thus, despite similar above- 
ground biomass yields in many cases, these species are likely to differentially affect soil C accumulation. 
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Introduction 


Recent attention focused on biomass crops to increase 
and diversify energy production and help mitigate 
greenhouse gas emissions has identified perennial 
grasses as potential dedicated energy feedstocks (Lemus 
& Lal, 2005; Carroll & Somerville, 2009; Davis et al., 
2010; Somerville et al., 2010; Don et al., 2012; Drewer 
et al., 2012). Although large plantings of perennial 
energy grasses could help to supplement biofuel pro- 
duction, the implications for other ecosystems services 
such as carbon (C) sequestration and nutrient cycling 
are not well understood. The decomposition of plant tis- 
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sues, particularly belowground pools, is an important 
process that affects soil microbial activity and C accu- 
mulation (Six et al., 2000; Jones & Donnelly, 2004; Rasse 
et al., 2005). As the main source of C inputs to the soil, 
the quantity and quality of root and rhizome biomass 
influences their decomposition rate and their residence 
time with subsequent impacts on C stored in the soil 
(Fontaine et al., 2007; Hattenschwiler & Jorgensen, 2010; 
Amougou et al., 2011; Knoll et al., 2012; Sun et al., 2013). 
Therefore, understanding the decomposition of below- 
ground biomass is critical in investigating C fluxes in 
terrestrial ecosystems. 

Root decomposition is greatly affected by climatic 
and edaphic conditions and is thus expected to be 
unique to each specific environment. In general, litter 
decomposition rates are positively correlated with mean 
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annual temperature and precipitation (Silver & Miya, 
2001; Zhang et al., 2008; Prescott, 2010). Similarly, satu- 
rated soil moisture conditions or extreme low moisture 
levels tend to limit decomposition of plant litter (von 
Haden & Dornbush, 2014; Lee etal., 2014). Roots 
decompose faster in clay loam than in sand and clay 
soils (Silver & Miya, 2001). 

Root chemical and morphological characteristics also 
contribute to their decomposition patterns and rates 
(Silver & Miya, 2001; Puttaso et al., 2011; Aulen et al., 
2012). Labile compounds in root tissues, such as simple 
sugars, organic acids, small-chain fatty acids, and pro- 
teins, can be easily taken up by microorganisms and 
metabolized, and their concentrations in the plant tissue 
tend to decline quickly as the decomposition process 
progresses (Rasse et al., 2005; Berg & McClaugherty, 
2008). Conversely, cellulose is a relatively recalcitrant 
cell wall component compared to water-soluble com- 
pounds and the decomposition of plant-derived C typi- 
cally decreases as cellulose concentration increases 
(Birouste ef al., 2012). However, research has also 
demonstrated no correlation between cellulose concen- 
tration and root decomposition rates (Aulen etal., 
2012). Lignin, suberin, cutin, and polyphenols are con- 
sidered recalcitrant components and generally retard 
root decomposition (Rasse et al., 2005; von Luetzow 
et al., 2006; Hattenschwiler & Jorgensen, 2010). In addi- 
tion, herbaceous species root morphological attributes 
such as specific root length (SRL) have been positively 
correlated with root decomposition rate (Aulen et al., 
2012). In a 4-year experiment with temperate tree spe- 
cies, root decomposition rates differed among root 
diameter classes (e.g., <0.5 and 0.5-2 mm) (Sun et al., 
2013). In contrast, Birouste et al. (2012) found no corre- 
lation between initial SRL or diameter and decomposi- 
tion rates. Consequently, there are conflicting reports in 
the literature on the effects of root chemistry and mor- 
phology on root decomposition rates, and the complex- 
ity is probably due to the interaction of root 
characteristics and environmental factors. 

In general, Poaceae roots decompose more slowly 
than those of herbaceous species from other taxonomic 
groups (Aulen et al., 2012; Birouste et al., 2012). The 
slow decomposition rates in Poaceae roots were often 
associated with their unique root morphological and 
chemical traits compared with other species, such as 
low N and P concentrations, high C/N, and high tissue 
density (Birouste et al., 2012; Siqueira da Silva et al., 
2015). Within the Poaceae family, species also differed 
in their root morphological and chemical traits and thus 
decomposition rates (Fornara & Tilman, 2008; Aulen 
et al., 2012). In addition, perennial warm-season grasses 
generally exhibit abundant root biomass and deep root- 
ing systems (Monti & Zatta, 2009; Somerville et al., 


2010). For instance, root dry weight in a 0- to 120-cm 
soil profile was 14 Mg ha! in giant reed (Arundo 
donax), 7.6 Mg ha! in giant miscanthus (Miscant- 
hus x giganteus), and 8.5 Mg ha™' in switchgrass (Pan- 
icum virgatum) in north Italy (Monti & Zatta, 2009). The 
large differences in root biomass production and root 
characteristics among perennial grass species are 
expected to result in differences in the amount and 
characteristics of C inputs to the soil. 

Despite relatively similar aboveground biomass pro- 
duction (Erickson et al., 2012; Knoll et al., 2012; Fedenko 
et al., 2013), our understanding regarding the potential 
impacts of root morphological characteristics and chem- 
ical composition of perennial bioenergy crops on below- 
ground C fluxes remains limited, particularly in the 
southeastern USA, where the warm temperatures asso- 
ciated with abundant rainfall are expected to create 
favorable conditions for root decomposition. Therefore, 
the objectives of this research were to characterize root 
morphological and chemical differences in six bioenergy 
grass species and to evaluate their impacts on root 
decomposition over a 12-month period under field con- 
ditions. We hypothesized that variability in species root 
chemistry and morphology would lead to differences in 
decomposition that will help to better understand C 
cycling in bioenergy grass cropping systems. 


Materials and methods 


Plant materials, growth conditions, and treatments 


Giant reed (Arundo donax), elephantgrass (Pennisetum pur- 
pureum ‘Merkeron’), energycane (Saccharum spp. ‘L79-1002’), 
sugarcane (Saccharum spp. ‘CP89-2143’), sweetcane (Saccharum 
arundinaceum ‘IK76-110’), and giant miscanthus (Miscant- 
hus x giganteus) were established from vegetative propagules 
from November 2008 to January 2009 at the University of Flor- 
ida Plant Science Research and Education Unit (29°24’N 
82°10'W) in Citra, Florida. The six species were arranged in a 
randomized complete block design with four replicates, giving 
24 plots in total. Plots consisted of six 6-m-long rows initially 
planted with 1-m row spacing (i.e, each plot was approxi- 
mately 36 m°). All plots were fertilized at a rate of 280 kg N 
ha~! yr~! split into applications of 90 kg N ha! in mid-April 
and 190 kg N ha™' in June. Plots also received ~ 70 and 
140 kg ha~’ of P,O; and K,O in each year, respectively, along 
with micronutrients. Limited irrigation was applied via over- 
head irrigation during establishment and with appearance of 
visible water stress (i.e., leaf rolling). Soil temperature at 10 cm 
was monitored during the experiment and is summarized in 
Fig. 1. 

For aboveground dry biomass yield, plots were harvested 
once annually in 2009 and 2010 (November) and twice annually 
(July and November) in 2011 and 2012. A 4-m section (4 m?) 
from one of the inner two middle rows in each plot was har- 
vested and weighed in the field. A representative subsample of 
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Fig. 1 Daily maximum (black line) and minimum (gray line) 
soil temperature (°C) at a 100-mm soil depth over the duration 
of the field root decomposition study. 


biomass was then collected and oven-dried at 50°C to a 
constant weight to determine dry matter concentration and 
dry biomass yield. 


Root morphology, chemistry, and decomposition 


To obtain a representative root biomass sample, at least four 
soil cores (11 cm diameter x 20 cm depth) were collected from 
each of the six perennial grass species plots after the above- 
ground biomass harvest in December 2011. Soil cores from the 
same plot were pooled into a composited sample and then 
taken back to the laboratory. Roots were separated from soil 
using a 2-mm sieve, gently washed with deionized water, and 
dried in the oven at 30 °C for 1-2 h to approximately 40% 
moisture. The roots for a given plot were then cut into 5-cm 
pieces, thoroughly mixed, and homogenized, and 3.45 + 0.05 g 
[mean (n = 92) sample weight +1 SE at 40% moisture] was put 
into 15 x 20 cm, 250-um mesh nylon litter bags (Castillo et al., 
2010). A subsample of cut roots was placed in the oven at 
50 °C and dried to a constant weight to determine the root 
moisture content and root dry weights for each bag. An addi- 
tional fresh root subsample from each plot for each species 
(n = 4) was scanned and analyzed with a digital image analysis 
system (WinRHIZO, Regent Instrument, Quebec, CA, USA) to 
determine root length, surface area, and volume. The subsam- 
ple was then placed in the oven at 50 °C and dried to a con- 
stant weight to determine root dry weight. Specific root length 
(SRL), specific root area (SRA), and specific root volume (SRV) 
were calculated as the ratio of root length, surface area, and 
volume divided by the corresponding root dry weight, respec- 


tively. 

Four root decomposition bags were buried horizontally at a 
7.5-cm depth in the soil in the row middles (50 cm away from 
each row) in their corresponding field plots on December 8, 


2011. For species without obvious plant rows after 2 years since 
planting, for example, giant reed, decomposition bags were 
buried arbitrarily in the center of the plots. One bag from each 
plot was collected at 1, 3, 6, and 12 months after installation. 
Root biomass remaining in the bags after decomposition was 
oven-dried at 50 °C to a constant weight and then was ashed at 
500 °C for at least 6 h to calculate the root biomass on an ash- 
free basis. The percentage of the initial root biomass remaining 
after decomposition (M, %) was calculated as follows: 


_ Rootmass,t = Ashmnass,t 
Rootmass,i = Ashmass,i 


t x 100, (1) 
where Rootmass,; and Rootmassy are the initial dry root biomass 
before decomposition and remaining at each harvest, respec- 
tively, and Ashmass,; and Ashmass, are the ash concentrations in 
the initial and remaining root biomass, respectively. 

To estimate the decomposition rate for each species, the pro- 
portion of the initial biomass remaining over time (t) was fit 
with a single-pool negative exponential model (Adair et al., 
2010; Birouste et al., 2012): 


M; = 100 x e*', (2) 


where K is the decomposition rate constant and is expressed in 
gkg! day! 

Root biomass from each plot before and after the 12-month 
field incubation was ground using a Restsch Mixer Mill 
(MM400, Verder Scientific, Inc., Newtown, PA, USA). A sub- 
sample of 30-50 mg of ground root was wrapped in tin cap- 
sules (9 x 10 mm, Costech Analytic Technology, Inc., Valencia, 
CA) and then the total C and N concentrations were measured 
by dry combustion using an elemental analyzer (FLASHEA 112 
Series, Thermo Fisher Scientific, Inc., Waltham, MA, USA). On 
another subsample, nonstructural extractives, structural carbo- 
hydrates, and lignin in root tissue before and after 12 months 
of decomposition were determined according to the procedures 
fully described in Fedenko et al. (2013). Briefly, to remove non- 
structural extractives, 0.5 g of ground root tissue was auto- 
claved with 100 ml of deionized water in 140-ml sealed 
pressure tubes (ACE Glass, Inc., Vineland, NJ, USA) at 121 °C 
and 103 kPa for 1 h. Autoclaved samples were then vacuum- 
filtered through coarse-porosity filter paper (> 25 um, What- 
man 113, GE Healthcare UK Limited, UK) to capture all struc- 
tural fiber. Captured structural biomass was then dried at 
50 °C to a constant weight for subsequent fiber carbohydrate 
and lignin analysis. A subsample of 0.3 g of structural fiber 
was incubated at 30 °C in 3 ml of 72% sulfuric acid for 1 h fol- 
lowed by digestion in 87 ml of 4% sulfuric acid (by adding 
84 ml deionized water) in an autoclave at 121 °C and 103 kPa 
for 1 h. Hydrolyzed samples were vacuum-filtered through a 
medium-porosity filtering crucible (Coors #60531, CoorsTex, 
Golden, CO, USA), and the filtrate was collected and analyzed 
for acid-soluble lignin using a UV-vis spectrophotometer (Stel- 
larNet, Inc., Tampa, FL, USA) at a wavelength of 240 nm. The 
filtered solids were dried at 105 °C to constant weight and then 
were ashed at 500 °C for at least 6 h. The residuals were 
weighed for ash concentration, and acid-insoluble lignin was 
calculated as the difference between dried filtered solids and 
ash. Total lignin (TL) was calculated as the sum of acid-soluble 
and acid-insoluble lignin. The remaining filtrate was adjusted 
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to pH to 5-7 with calcium carbonate and filtered through a 
0.22-um syringe filter (Fisher Scientific, Pittsburgh, PA, USA) 
for fiber carbohydrate analysis using high-performance liquid 
chromatography (HPLC). The filtrate samples were analyzed 
by HPLC (Perkin-Elmer Flexar system, Waltham, MA) with an 
Aminex HPX-87H column (Bio-Rad, Hercules, CA) maintained 
at 50 °C with HPLC-grade 4 mM sulfuric acid as the mobile 
phase at 0.4 ml min™ with a 10-4] injection and 40-min run 
time. Concentrations of total fiber glucose (TFG), total fiber 
xylose (TFX), total fiber arabinose (TFA), total fiber carbohy- 
drates (TFC) (the sum of TFG, TFX, AND TFA), and TL were 
all expressed on an ash-free basis (mg g~' DM). 

The remaining chemical components after the 12-month field 
incubation, including total lignin (TLremaining), total fiber glu- 
cose (TFG,emaining), total fiber xylose (TFX;emaining), total fiber 
arabinose (TFA;emaining), total fiber carbohydrate (TFCyemaining), 
nitrogen (Nremaining), and carbon (Cyemaining), were calculated 
using the following equation (Fioretto et al., 2005): 


cf X Rootmass, 
= ea ae x 100, (3) 
where Rc is the remaining chemical component after the 
12-month incubation in the field. C; and C; are initial and final 
concentrations of each chemical component, respectively. 
Rootmass,s is the final weight of root biomass after the 12-month 
incubation. 


Data analysis 


Root chemical (e.g., concentrations of TL, TFG, TFX, TFA, TFC, 
and total C and N) and morphological (e.g., SRL, SRA, and 
SRV) characteristics prior to decomposition were analyzed 
using the generalized linear mixed model (glimmix) procedure 
of SAS (ver. 9.3, SAS institute, Cary, NC, USA) with species as 
a fixed effect and block as a random effect. Remaining chemical 
components and K were also analyzed using the glimmix pro- 
cedure. Pairwise comparisons were made using the Ismeans 
statement with the Tukey’s method with a significance level of 
P < 0.05. Relationships between K and root traits were ana- 
lyzed with Pearson’s correlation, principal component analysis 
(PCA), and factor analysis using SAS procedures of corr, prin- 
qual, and factor, respectively. 


Results 


Decomposition rate 


The decomposition rate constant, K, of elephantgrass 
was higher at 3.64 g kg ' day ' than all other species 
except for sweetcane (Table 1). This was associated with 
23% of root biomass remaining after the 12-month incu- 
bation (Fig. 2), which was also among the lowest of all 
species in this study. Sweetcane decomposed at a rate 
constant of 2.77 g ke! day ', which was significantly 
higher than sugarcane and energycane, but did not dif- 
fer from elephantgrass, giant reed, or giant miscanthus. 
After the 12-month decomposition study, 40% of sweet- 


cane root biomass remained, which was lower than that 
in giant reed and energycane (data not shown). Energy- 
cane, sugarcane, giant reed, and giant miscanthus did 
not differ in K, averaging 1.72 g ke! day '. 


Root morphology 


Specific root length varied almost three-fold among spe- 
cies, ranging from 9.0 for giant reed to 25.3 m gs! for 
elephantgrass (Table 1). Elephantgrass roots possessed 
the highest SRL, whereas giant reed, sugarcane, sweet- 
cane, and energycane were among the lowest in SRL. In 
contrast, there was less than two-fold variation in SRA, 
which ranged from 199 to 337 cm? g`. Elephantgrass 
and giant miscanthus roots had the highest SRA, while 
giant reed, sweetcane, and sugarcane were among the 
lowest in SRA. Specific root volume ranged from 3.06 to 
4.76 cm? g`" across all species. Giant miscanthus, ener- 
gycane, and sugarcane roots were among the highest in 
SRV, while elephantgrass, sweetcane, and giant reed 
were among the lowest for SRV. 


Root chemical composition 


Elephantgrass roots exhibited a lower TE concentration 
prior to the 12-month incubation than giant miscanthus, 
which was not different from other species (Table 2). It 
also had lower TL concentration than giant reed. Root 
concentrations of TFG and TFX in elephantgrass roots 
were higher than almost all other species, in contrast to 
the relatively lower concentrations in sugarcane. Root 
concentrations of TFA were relatively lower in giant 
reed, sweetcane, and elephantgrass. Overall, root TFC 
concentration of elephantgrass was higher than all other 
species except sweetcane. Root total C varied little 


Table 1 Rate constant (K) for decomposition, specific root 
length (SRL), specific root surface area (SRA), and specific root 
volume (SRV) of perennial grass fine roots before decomposi- 
tion 


K SRL SRA SRV 
Species gkg day! meg! cmg! cme! 
Giant reed 1.80 BC* 90C 204 C 3.59 BC 
Sugarcane 1.62 C 10.5C 231 BC 3.97 AB 
Sweetcane 2.77 AB 10.1 C 199 C 3.06 C 
Energycane 1.48 C 13.1 BC 259 B 4.28 AB 
Elephantgrass 3.64 A 253A 311A 3.09 C 
Giant 1.99 BC 18.9 B 337 A 4.76 A 
miscanthus 
P-value <0.001 <0.001  <0.001 <0.001 


*Means (n = 4) followed by different letters within a column 
differ significantly (P < 0.05). 
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Fig. 2 Root decomposition (% remaining root mass) over the 12-month decomposition period. Species include giant reed (GR), sug- 
arcane (SC), sweetcane (SW), energycane (EC), elephantgrass (EG), and giant miscanthus (GM). Circles represent treatment mean 
(n = 4) fraction of root biomass on an ash-free basis remaining at each removal date, and the error bars represent +1 SD. Solid lines 
represent the fitted lines from the single-pool negative exponential model, M, = 100 x e‘*”, for each of the species, where K is the 


decomposition rate in g kg~' day! and t is time in day. 


Table 2 Concentrations of total extractives (TE), total lignin (TL), total fiber glucose (TFG), total fiber xylose (TFX), total fiber arabi- 
nose (TFA), total fiber carbohydrate (TFC), total C, total N, C:N, and TL:N ratios of root tissue on an ash-free dry matter basis prior to 


decomposition 

Species TE TL TFG TFX TFA TFC Total C Total N C:N TL:N 
mg g ' DM 

Giant reed 152 AB* 327A 269BC 168 BC 19.4 C 455 BC 408B 61C 67.2 AB 53.8 A 

Sugarcane 146 AB 316 AB 248C 162 C 32.3 B 442 C 439 AB 8.0 AB 54.9 BC 39.6 BC 

Sweetcane 144 AB 299 AB 301B 195 AB 23.6C 519 AB 428AB 56C 76.6 A 53.6 A 

Energycane 155 AB 301 AB 267BC  176BC 387A 481 BC 432AB 6.8 BC 63.7 AB 44.3 AB 

Elephantgrass 126 B 293 B 342 A 207 A 20.6 C 570 A 456 A 7.8 AB 58.7 BC 37.7 BC 

Giant miscanthus 174 A 291 B 253C 172 BC 35.2 AB 459BC 416AB 88A 48.6 C 34.0 C 
0.058 0.003 <0.001 0.001 0.001 <0.001 0.048 <0.001 <0.001 <0.001 


P-value 


*Means (n = 4) followed by different letters within a column differ significantly (P < 0.05). 


among species (408-456 mg g ' DM); however, total C 
of elephantgrass was higher than that of giant reed. 
Root total N concentration ranged from 5.6 mg g | in 
sweetcane to 8.8 mg gt in giant miscanthus. Giant mis- 
canthus, sugarcane, and elephantgrass were among the 
species with the highest root N concentrations. Root C: 
N ratio ranged from about 49 in giant miscanthus to 
about 77 in sweetcane. Sweetcane, energycane, and 
giant reed were among the species with the highest root 
C:N ratios. The high TL:N ratio in giant reed resulted 


from its high TL concentration and low total N concen- 
tration. However, the relatively rapidly decomposing 
elephantgrass did not show lower C:N or TL:N ratio 
than most other species. 

Root chemical traits were analyzed again after their 
12-month incubations and expressed as a fraction of 
their initial amounts to better understand relevant dif- 
ferences in decomposition (Table 3). Elephantgrass roots 
exhibited relatively lower remaining TFG, TFX, TFA, 
TFC, TL, total C, and N (Table 3), which was consistent 
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with its higher K (Table 1). In contrast, roots of giant 
reed, sugarcane, and energycane showed higher remain- 
ing amount of all chemical characteristics (Table 3), 
which was consistent with their lower K (Table 1). 

Among all the six species, initial SRV was positively 
correlated with remaining TL, TFC, C, and N (Table 4). 
Remaining TL, TFC, C, and N were also positively cor- 
related with each other. 


Relationships between root decomposition rate and root 
traits 


K was positively correlated with a number of root 
chemical traits prior to decomposition, including TFG, 
TFX, TFC, and total C, but it was negatively correlated 
with TL and TFA concentrations and SRV (Table 5). 
Root TL was negatively correlated with total N concen- 
tration, SRL, and SRA. However, TFC concentration 
was positively correlated with SRL and negatively cor- 
related with SRV. 

The first two axes of the PCA performed with 12 root 
traits and K accounted for 78% of the variance (Fig. 3). 
The first PCA axis (Component 1) accounted for 42% of 
the variance and was defined by root chemical and mor- 
phological traits. The second PCA axis (Component 2) 
accounted for 36% of the variance and was defined by 


K. Concentrations of TFG, TFX, TFC, and total C were 
grouped together with K, indicating a high correlation 
among these traits. Additionally, SRL was slightly posi- 
tively correlated with K. However, SRV, TFA, and TL in 
the root tissue predecomposition were negatively 
related to K. Root total N concentration, C:N, TL:N, and 
SRA prior to decomposition were independent of K, as 
indicated by the near-90° angle among the directional 
vectors. 


Discussion 


Root chemistry plays a dominant role in controlling pat- 
terns of decomposition rates at a global scale (Silver & 
Miya, 2001). Previous studies have commonly focused 
on plant tissue C:N ratios, which were often negatively 
correlated with decomposition rate as described in 
Zhang et al. (2008). However, decomposition is more 
complex, and rates of decomposition are not necessarily 
related to C:N ratio, but they can be influenced by a 
number of factors such as environmental conditions, 
decomposer composition, other root chemical compo- 
nents, and physical structure (Table 5) (Johnson et al., 
2007; Birouste et al., 2012; Smith et al., 2014). 

Although there are a number of factors that have 
been shown to influence decomposition, there is a grow- 


Table 3 Amount (% of initial) of total fiber glucose (TFGyemaining), total fiber xylose (TFX;emaining), total fiber arabinose (TFA; emaining), 
total fiber carbohydrate (TFCyemaining), total lignin (TLyemaining), C (Cremaining), and N (Nyemaining) remaining on an ash-free dry matter 
basis after a 12-month decomposition period in the field for the six perennial grass species 


Species TFG, emaining TEXyemaining TFA yemaining TFCyemaining TLyemaining Cyemaining Nyemaining 
Giant reed 50 AB* 48 AB 58 A 50 AB 74 AB 58 A 76 AB 
Sugarcane 53 AB 51 AB 60 A 53 AB 70 AB 55 AB 83 AB 
Sweetcane 28 BC 28 BC 34 B 28 BC 47 BC 33 BC 57 BC 
Energycane 61 A 54 A 66 A 59 A 76 A 63 A 91 A 
Elephantgrass 16C 16C 17 B 16 C 26 C 17C 28 C 
Giant miscanthus 34 BC 34 BC 34 B 34 BC 54 A-C 42 AB 60 B 
P-value 0.0004 0.0007 <0.0001 0.0004 0.0010 0.0001 0.0001 


*Means (n = 4) followed by different letters within a column differ significantly (P < 0.05). 


Table 4 Pearson’s correlation coefficients between remaining total lignin (TLyemaining), fiber carbohydrate (TFC,emaining), C (Cremaining) 
and N (Nyemaining) after 12-month decomposition, and specific root length (SRL), surface area (SRA), and volume (SRV) before decompo- 


sition of all perennial grass species 


TLyemaining TFC, emaining Cyemaining Nyemaining SRL SRA 
TFC, emaining 0.84** ai 
Cremaining CIZER 0.837** 
Nyemaining 0.93** 0.82*** 0.94*** 
SRL —0.34 —0.40 —0.26 —0.41* 
SRA —0.05 —0.14 —0.02 —0.15 0.87*** 
SRV 0.52** 0.47* 0.52** 0.44* —0.00 0.48** 


+P < 0.01, 0.05, and 0.1 represented as ***, **, and *, respectively. 
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TL:N SRL SRA 


C:N 


Total N 
—0.94%*** 
_0.95%* 


Total C 
0.13 
0.17 


TFC 
0.46** 

—0.16 
0.30 


—0.31 
—0.02 
0.38* 
—0.38* 
—0.40* 
0.00 
0.36 


TFA 


—0.29 
0.99"** 
0.40* 

—0.17 
0.31 


TFX 


—0.51** 
0.97*** 
0.45** 

—0.21 
0.34 
0.13 


TFG 


—0.17 
—0.22 
—0.19 
—0.23 
—0.15 
—0.39* 
0.32 


TL 


0:757** 
0.63*** 
_0.64*** 
0.67*** 
0.40* 
—0.05 
0.10 


Table 5 Pearson’s correlation coefficients between root decomposition rate constant (K) and root morphological and chemical traits among all perennial grass species. Root traits 
—0.45**+ 


prior to decomposition include concentrations of total lignin (TL), total fiber glucose (TFG), total fiber xylose (TFX),total fiber arabinose (TFA), total fiber carbohydrate (TFC), total 


C and N in root tissue, ratios of TL:N and C:N, and morphological traits of specific root length (SRL), specific root surface area (SRA), and specific root volume (SRV) 


Total C 
Total N 


TFA 
TFC 


Z 
Ü 


TL 
TFG 
TF. 


—0.18 0.90*** 
—0.51** 


0.06 


0.09 


0.60*** 
—0.48** 
—0.44* 


—0.15 
0.38 
0.08 


Z 
4 


_0.61*** 
_0.69%** 
—0.35 


0:52** 


0.13 
—0.06 
—0.36 


0.49** 
0.13 
—0.66*** 


0.47** 
0.12 
—0.66*** 


0.45** 
0.05 


SRL 


0.87*** 
—0.00 


_0.67*** 
—0.52** 


0.63*** 
0.41* 


5 


0.48** 


0.03 


—0.58** 


SRV 
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+P < 0.01, 0.05, and 0.1 are represented as ***, **, and *, respectively, which is the same in the following tables for correlation analysis. 


SRL 


SRA 
Total N 


NE 


Component 2 (36.02%) 


| C | | 1 | | | | 
-2 -1 0 1 2 3 
Component 1 (41.85%) 


Fig. 3 Principal component analysis for root traits and decom- 
position rate constant (K). Root traits prior to decomposition 
include specific root length (SRL), specific root surface area 
(SRA), specific root volume (SRV), and concentrations of total 
fiber glucose (TFG), total fiber xylose (TFX), total fiber arabi- 
nose (TFA), total lignin (TL), total C and N in root tissue, and 
ratios of TL:N and C:N. Species include giant reed (GR), sugar- 
cane (SC), sweetcane (SW), energycane (EC), elephantgrass 
(EG), and giant miscanthus (GM). 


ing body of literature indicating that tissue carbohy- 
drate concentration and composition are key determi- 
nants of root decomposition rates. For example, high 
fiber carbohydrate concentrations were associated with 
rapid root decomposition among tree and herbaceous 
species (Silver & Miya, 2001; Aulen et al., 2012). Results 
from the present study confirm these findings, as root K 
was positively correlated with TFG, TFX, and TFC prior 
to decomposition among botanically similar (i.e., Poa- 
ceae family) perennial grass species (Table 5 and Fig. 3). 
The positive correlation between TFC and K, and in par- 
ticular, TFG and K, could be explained by energy domi- 
nance. Decomposer activity during bio-decomposition is 
mainly controlled by the energy that can be supplied by 
substrates contained in the litter (Fioretto et al., 2005; 
Hattenschwiler & Jorgensen, 2010). Compared with lig- 
nin, cellulose is an energy-rich compound. Thus, the 
faster decomposition of elephantgrass roots in the pre- 
sent study could be because of their relatively higher 
cellulose concentration and _ cellulose:lignin ratio 
compared with the other species. 

On the other hand, the energy supplied by more 
recalcitrant compounds is generally low and the energy 
required to break them down is high (Fioretto et al., 
2005; Hattenschwiler & Jorgensen, 2010; Sun et al., 
2013). The energy is therefore not enough to support the 
level of microbial activity that is seen with less recalci- 
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trant compounds, as indicated by the negative 
correlation between K and initial TL in the present 
study (Table 5). Thus, the decomposition of recalcitrant 
compounds was relatively slow, as indicated by the 
generally high remaining TL after the 12-month field 
incubation, with the exception of elephantgrass 
(Table 3). In the present study, elephantgrass possessed 
the highest initial TFC concentration among all the six 
perennial grass species, but its initial TL was similar to 
other species except for giant reed (Table 2). Further- 
more, consistent with the root decomposition study in 
the field, remaining root biomass was also the lowest in 
elephantgrass among the six species in another decom- 
position study performed in pots without plants inside 
a greenhouse over a 12-month period (data not shown). 
Root K of elephantgrass from the present study is very 
close to that from another study with K values ranging 
from 2.65 to 3.18 g kg ' day ' (Siqueira da Silva et al., 
2015). Taken together, rapid decomposition in elephant- 
grass root tissue, including TL, in the present study 
could be explained, at least in part, by a relatively high 
TFC and TFG to TL ratio. The TFC and TFG provided 
sufficient energy to decompose more recalcitrant com- 
pounds similar to the priming effect seen with root exu- 
dates and fine root turnover (van der Krift et al., 2002; 
Fioretto et al., 2005; Talbot & Treseder, 2012). 

As a result, the decomposition of labile and recalci- 
trant compounds was correlated in the study, as indi- 
cated by the positive correlation between TFCyemaining 
and TLyemaining (Table 3). In the present study, elephant- 
grass possessed the highest initial TFC concentration 
among all the grasses, but its initial TL was similar to 
other species except for giant reed (Table 2). Its high ini- 
tial TFC concentration probably stimulated or primed 
lignin decomposition, which resulted in relatively low 
amounts of both TFCyemaining (16%) and TLyemaining 
(26%) after the 12-month field incubation. On the other 
hand, giant reed had at least 50% of both TFC and TL 
remaining after the incubation. Consequently, a rela- 
tively high TL and/or high TL:TFC reduced K by slow- 
ing the loss of hemicellulose and cellulose as well as 
lignin, implying that lignin may have protected the cell 
wall from degradation (Fioretto et al., 2005; Talbot & 
Treseder, 2012). 

Besides the physical protection of more easily decom- 
posable components (i.e., cellulose and hemicellulose), 
the recalcitrance of lignin is related to its chemical char- 
acteristics, such as molecule size, polarity, three-dimen- 
sional structure, and functional groups (aromatic ring 
structures) (von Luetzow et al., 2006; Puttaso et al., 
2011; Gul & Whalen, 2013). Once plant tissues start to 
decompose, lignin begins to incorporate N, and conden- 
sation reactions take place (Berg & McClaugherty, 
2008). These chemical transformations could cause 


changes in structures that are resistant to degradation 
by soil microbes and also act as barriers limiting their 
access to the more labile compounds. High N concentra- 
tions also suppress the formation of ligninase (Berg & 
McClaugherty, 2008) and may help to explain why a 
positive correlation was found between TLremaining and 
N, which was consistent with the previous findings, 
especially in lignin-rich plant tissues (Hobbie, 2000; Per- 
akis et al., 2012). Also, because of the protection of lig- 
nin, degradation of TFC could have been retarded, 
leading to the positive correlation between remaining 
TFC and N (Table 4). 

Root morphological traits represent the economics of 
root investment: carbon input for root growth vs. the 
capacity of resource acquisition (Donovan et al., 2014; 
Reich, 2014). Their correlations with root decomposition 
rate reflect their potential in C and nutrient cycling 
across ecosystems (Donovan et al., 2014). In the present 
study, morphological traits associated with fast root 
decomposition were high SRL and, even more so, low 
SRV roots (Fig. 3; Table 5). These findings were in 
agreement with other studies that have shown SRL to 
be positively correlated with decomposition rate 
because a high SRL has the potential to facilitate decom- 
position through maximizing root surface area and 
exposure for bio-decomposition (Aulen et al., 2012; Bir- 
ouste et al., 2012; Donovan et al., 2014; Smith et al., 
2014). In addition, the present study identified correla- 
tions between root morphological and chemical traits, 
such as the positive correlation between TLyemaining and 
SRV (Table 4), and the negative correlations between 
TL, SRL, and SRA (Table 5). Root diameter has also 
been shown to affect root decomposition rates among 
switchgrass cultivars (de Graaff etal., 2013), which 
could be related to different concentrations of soluble 
carbohydrates and lignin in different root diameter size 
classes (Fan & Guo, 2010). However, interactive effects 
of root morphological and chemical traits on root 
decomposition have not been investigated thoroughly 
in the previous research. Varied morphological and 
chemical traits of plant species are determined by both 
genetic and environmental components. The interaction 
of root morphological and chemical characteristics, thus, 
might play an important role in root decomposition 
under different environments. For instance, after similar 
time periods of decomposition in the field, the remain- 
ing root biomass of the same species can differ greatly 
under varied environments (Harmon et al., 2009). 

Beyond the intrinsic characteristics of roots, environ- 
mental factors were also likely to influence root decom- 
position. For instance, soil temperature has been shown 
to be closely correlated with root decomposition rates, 
as decomposition rates increased up to five-fold with a 
temperature increase from 20 to 30 °C (Solly et al., 2014; 
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Stewart et al., 2015). In the present study, root mass loss 
was accelerated at 185 days after decomposition (Fig. 2), 
which was consistent with the high soil temperature in 
June (Fig. 1). Additionally, species or environmental 
factors (e.g., soil temperature, moisture, and fertility) 
and their interactions could have contributed to differ- 
ences in microbial community diversity and/or func- 
tional activity (van der Heijden et al., 2008; Brzostek 
et al., 2015) that could have contributed, at least in part, 
to the observed differences in root decomposition 
among the species as well. 

Overall, perennial grass species, many of which pos- 
sess similar aboveground annual biomass production, 
differed substantially in root morphological and chemi- 
cal traits, and these traits were closely correlated with 
their decomposition rates. In the context of C flux across 
ecosystems, root biomass (i.e., quantity and quality) and 
its decomposition rate are among the most important fac- 
tors governing soil C dynamics (Aerts et al., 2003; von 
Luetzow et al., 2006). Although the processes controlling 
belowground C translocation and allocation of root- 
derived C into soil are complex, more organic matter 
input from roots combined with greater microbial activ- 
ity are presumed to lead to more C transfer and storage 
in the soil (Kuzyakov & Domanski, 2000; Puttaso et al., 
2011). However, the quality of root-derived C inputs can 
also play a dominant role in the transfer of root biomass 
to soil organic C. Despite a lack of consensus in the litera- 
ture, most recent studies suggest that high-lignin materi- 
als cannot be used efficiently by the soil microbial 
community, and high-lignin materials thus contribute 
less to soil organic C than materials with low lignin con- 
centrations (Hancock et al., 2007; Cotrufo et al., 2013; Ste- 
wart etal., 2015). Additionally, introducing readily 
available organic matter into soils initially brings up 
priming effects, which causes soil organic C loss through 
promoted microbial metabolisms (Brzostek et al., 2015; 
Stewart et al., 2015). However, the negative effects of 
root-derived C inputs on soil C stocks are generally tran- 
sient and more pronounced in the rhizosphere (Haichar 
et al., 2014). As evidenced in previous studies, the extent 
that root-derived C affects soil organic C in bioenergy 
production systems is expected to vary considerably 
depending on soil, plant, and environmental factors 
(Bandaru et al., 2013; Bonin & Lal, 2014). Thus, changes 
in soil organic C can be expected to differ between bioen- 
ergy grasses. For example, after conversion of native 
land to miscanthus, soil organic C was building up 
slowly over time, whereas the conversion to sugarcane 
caused a large initial loss of soil organic C in the top soil, 
and the loss could last for a few decades before soil 
organic carbon rebuilds (Anderson-Teixeira et al., 2009). 
Such interspecies differences in soil organic C could be 
related to the characteristics of root chemistry and 


morphology indicated in the present study. For instance, 
giant miscanthus and giant reed showed slow root 
decomposition rates, implying a slow increase in soil 
organic C. In contrast, elephantgrass roots exhibited all 
the root traits favorable for fast decomposition: high 
TFG, high TFC, high SRL, and a low SRV, and elephant- 
grass could thus be expected to contribute to soil C accu- 
mulation in a rapid way. 
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